AC driven thermal ratchets. 
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Abstract : 

We consider the motion of a overdamped Brownian particle in periodic asymmetric potential with space dependent friction 
coefficient. In the presence of external time periodic forcing, the system shows multiple current reversals on varying the 
amplitude of the external forcing and the temperature of the thermal bath. In the adiabatic regime we find a single reversal of 
current as a function of noise strength which can only be accounted due to the presence of space dependent friction coefficient. 
For very large forcing term, the current does not go to zero, instead it asymptotically tends to a limiting value depending on 
the phase shift between the potential and the friction. This fact plays an important role in obtaining multiple current reversals. 



1. Introduction 

Rectification of noise leading to unidirectional motion 
in ratchet systems has been an active field of research 
over the last decade. In these systems directed Brownian 
motion of particles is induced by nonequilibrium noise 
in the absence of any net macroscopic forces and po- 
\ tential gradients. Several models have been proposed to 

■ explain this transport mechanism under various nonequi- 
librium situation, like (a) Flashing ratchets [1,2], wherein 
the particles experience a fluctuating energy profile, (b) 
Rocking ratchets [3,4], where the Brownian particles expe- 

, riences a spatially uniform, time periodic uniform force, 
F(t + t) = F(t), (c) Diffusion ratchets[s], wherein par- 
| tides are driven by a time periodic diffusion coefficient 

■ D(t) = D(t + t) and (d) Correlation ratchets [6], where 
particles are driven by spatially uniform but temporally 
correlated force ((t) of zero average. In these systems to 

' get a unidirectional current either spatially asymmetric 
, periodic potentials or time asymmetric external forces are 

■ necessary. The recent burst of work on ratchet systems 
is motivated in part to explain the unidirectional trans- 
port of molecular motors in biological systems and a new 
prospect of novel techniques for separation or segrega- 
tion of particles of macro-meter size. This technique for 
particle separation is based on the property of current re- 
versal in ratchet by varying the strength of thermal noise, 
amplitude of the forcing, size of the particle or any other 
relevant variable in the problem. 

In the area of Brownian rectifiers or ratchets the study 



of current reversal has become a subject by itself. Bar- 
tussek et. al.[i] showed the occurrence of current rever- 
sals in rocked thermal ratchet as a function of amplitude 
of rocking force as well as the temperature of thermal 
bath. Multiple current reversals have also been shown in 
the deterministic limit of these ratchets when the iner- 
tial term is taken into account. However, these multiple 
current reversals in inertial ratchets are not robust in the 
presence of noise. Beside rocking ratchets, current rever- 
sals have also been observed in flashing ratchets[7,2]. We 
in this work report that that multiple reversals can be 
achieved even in rocked overdamped ratchet in the pres- 
ence of space dependent mobility, as a function of the 
noise strength and amplitude of the rocking force. In 
the deterministic overdamped case we get current rever- 
sal as a function of the amplitude of rocking force. It 
is to be noted that systems with space dependent fric- 
tion are not uncommon. Brownian motion in confined 
geometries show space dependent friction [s]. It is be- 
lieved that molecular motor proteins move close along 
the periodic structure of microtubules and will therefore 
experience a position dependent mobility[9]. Frictional 
inhomogeneities are common in super lattice structures 
and Josephson junctions etc. 



2. Model 

We consider an overdamped Brownian particle mov- 
ing in a asymmetric potential V(x) with space depen- 
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dent friction coefficient 77 (x) under the influence of ex- 
ternal force field -F(t) at temperature T. Through- 
out our analysis we take the ratchet potential V(x) = 
— ^:(sin(27ra;) + j sin(47ra;)). Here [i is the asymmetry 
parameter with values taken in the range < n < 1, 
friction coefficient r)(x) = 770(1 — A sin(27rx + </>)), |A| < 1. 
<t> determines the relative phase shift between friction co- 
efficient and potential . The forcing term is taken to 
be F(t) = Asin(w£ + 6), (ui = — , where r is period of 
force) . Without any loss of generality 9 is taken to be 
zero. The Langevin equation for this system in the over- 
damped limit is given by[ii,i2] 



4e-03 




m, (i) 



(V'(x)-F(t)) 
rj(x) 



where £(t) is Gaussian thermal noise with correlation 
< >= 28 (t - t') The Fokker Planck Equation 

( FPE)[n] corresponding to Eqn. (1) is given by 



dP(x,t) _ d 1 



dt 



dx rj{x) 



[k B T-^- + (V'(x)-F(t))}P(x,t). 



(2) 



Here P(x,t) is the probability density at position x 
at time t. Equation (2) can be recasted to a continuity 
equation — 7^- 
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J(x,t) = - — [(V'(x) - F(t)) + k B Tj-]P(x,t), (3) 

is the probability current density. Since the potential 
and the driving force have spatial and temporal peri- 
odicity respectively, therefore J(x,t) = J(x + l,t + r), 
[443]. The net current j in the system is given by 
j = lim^oo i J t * +T dt Jq J(x, t)dx . It should be noted 
that for symmetric potential and A = 0, j = 0. Rec- 
tification of current is possible if the potential is either 
asymmetric or A ^ with <fi ^ 0, ir [14]. It is impor- 
tant to note that in these systems rectification is due to 
a symmetry breaking as the potential is asymmetric and 
A = 0. For the case when the potential is symmetric and 
A =/= 0, the symmetry breaking results from the dynam- 
ics of the system. It is the space dependent mobility that 
breaks the symmetry in this noncquilibrium problem and 
the magnitude of the symmetry breaking is related to the 
phase shift between the friction and potential profile, j is 
independent of the initial phase 9 of the driving force. We 
solve the FPE numerically by method of finite difference 
and calculate the current j. All the physical quantities 
such as j, T, A, u> are in dimensionless units [13-15]. 



3. Results and Discussions 

In the Fig. (1A), the average current j is plotted as a 
function of temperature T for different values of u>. Here 
the asymmetry parameter n = 1.0 and A = 0.0. 
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Figure 1. Mean current j vs temperature T for (a) 
A = 0,/Lt = 1.0 and (b) A = 0.1, fx = and (f> = 0.2?r. Note 
there is no current reversal when the potential is symmetric. 



In this case the current reverses its sign (only once) for 
frequencies sufficiently large as shown in the u> = 4.0, u = 
5.0 case. In the absence of asymmetric potential and 
presence of space dependent friction ( A = 1.0), there is 
no current reversal irrespective of w and phase shift (f> 
as shown in Fig (IB). Hence asymmetric potential is es- 
sential for current reversal. However the direction of the 
current depends on the phase lag <f>. Separately in both 
these cases absolute value of current exhibits a maxima as 
a function of T, reminiscent of stochastic resonance phe- 
nomena. In a purely asymmetric case (A = 0) current 
vanishes rapidly when T exceeds the temperature associ- 
ated with the barrier height. Whereas, in the symmetric 
case due to space dependent friction absolute values of 
currents are significantly higher and decay slowly to zero 
in the large temperature regime. Naturally in the pres- 
ence of both asymmetry and space dependent friction for 
the case under study the low temperature regime is dom- 
inated by the effect of asymmetry while the high temper- 
ature regime is dominated by space dependent friction. 
From this , one can qualitatively explain the current re- 
versals from positive to negative side function of 
temperature even in the adiabatic limit. Hence in the 
presence of both space dependent friction and asymmet- 
ric potential we can have current reversal with T even in 
the adiabatic regime for a suitablily chosen value of <j) as 
shown in Fig (2A), where we have plotted j vs T with 
/1 = 1.0, A = 1.0 and <j> = 0.2-7T for different values of u>. 
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Figure 2. The mean current j vs temperature T for 
4> = 0.2-7T, A = 0.5 and A = 0.1. The driving frequencies 
are u = 3.0, 4.0 and 5.0. The right hand side figure shows 
current j vs T for (f> = 1.2n. 
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Figure 3. The mean current j with amplitude A of the 
forcing term for <j> = 0.88tt, T = 0.05 and A = 0.1 with 
lo — 3.0, 4.0, 5.0. The inset shows the reversal of determin- 
istic current vs the amplitude of the forcing. 



In the presence of finite frequency drive there are twice 
current reversals as shown in the figure. This phenom- 
ena of twice current reversal with temperature T is the 
foremost feature of our system, previously unseen in any 
overdamped system. When the phase difference <j> is such 
that the current due to space dependent friction alone is 
in the same direction as that of current due to potential 
asymmetry only, then we do not have current reversal in 
the adiabatic case as shown in Fig (2B) where <fi = 1.2tt, 
though a single current reversal due to finite frequency 
drive may be present. In all the cases studied so far 
we have observed that the current reversal don not take 
place above a critical frequency lo c of driving, which in 
turn depends on <f> and other parameters in the problem. 



In Fig. 3, the plot of j versus A is shown for differ- 
ent values of u, keeping A, <j) an d T fixed at 0.1, 0.887T 
and 0.05 respectively. For u> = 4.0 curve, we can see as 
many as four current reversals. For very large value of 
A, the current asymptotically goes to a constant value 
depending on the value of <j>, as was previously shown for 
the adiabatic case[i4]. The inset in Fig. 3 shows current 
reversal even for the deterministic case also. In addition 
it exhibits an interesting phenomena of current quanti- 
zation and phase locking. These features are washed out 
with increasing noise strength. 
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Multiple current reversals can also be seen when the 
amplitude ( A) of the forcing term is varied in a suitable 
parameter regime of our system. 



Figure 4. Mean current j vs temperature T for (a) 
A = 0, /u, = 1.0 and (b) A = 0.1, /i = and <f> = 0.2tt. Note 
there is no current reversal when the potential is symmetric. 
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In Fig (4) we have plotted j versus A for <fr = 1-2tt. 
There is no current reversal in the adiabatic regime for 
both deterministic as well as finite temperature case. 
The observation of multiple current reversals can be at- 
tributed to a cooperative interplay between the spatial 
asymmetry of the potential, the friction inhomogeneity 
and the finite frequency drive. Depending on the system 
parameters we may have multiple current reversal or no 
current reversal at all ( see Fig.( 3 and 4)). 

In conclusion, we have studied the transport properties 
of overdamped Brownian particles moving in an asym- 
metric potential with space dependent friction coefficient 
and rocked by periodic force. We observe multiple cur- 
rent reversal with both temperature T and A the am- 
plitude of the rocking force in the presence of finite fre- 
quency driving. Current reversal also occurs in the de- 
terministic adiabatic regime for suitable values of phase 
difference cf>. All the above results can be understood 
qualitatively. The space dependent friction plays an ac- 
tive role in these systems in determining several counter 
intuitive phenomena[i4-i6] observed in the transport pro- 
cesses. 
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